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Abstract 

We analyse a simple model of the heat transfer to and from a small satellite or- 
biting round a solar system planet. Our approach considers the satellite isothermal, 
with external heat input from the environment and from internal energy dissipation, 
and output to the environment as black-body radiation. The resulting nonlinear or- 
dinary differential equation for the satellite's temperature is analysed by qualitative, 
perturbation and numerical methods, which show that the temperature approaches 
a periodic pattern (attracting limit cycle). This approach can occur in two ways, 
according to the values of the parameters: (i) a slow decay towards the limit cycle 
over a time longer than the period, or (ii) a fast decay towards the limit cycle over a 
time shorter than the period. In the first case, an exactly soluble average equation is 
valid. We discuss the consequences of our model for the thermal stability of satellites. 

1 Introduction 

The design process of any spacecraft necessarily includes a thermal analysis, performed to 
guarantee that both the bus and the payload are kept along the whole mission within the 
appropriate temperature range [I]. The thermal environment of a satellite in orbit is very 
harsh, dominated by vacuum conditions which make difficult to find a way of cooling down 
the components of the spacecraft. Every single element of the satellite has to be thermally 
analysed to verify that its requirements are fulfilled. To do it, it is necessary a rather 
detailed design of the spacecraft, as it implies to calculate thermal loads and couplings 
between elements. This is usually carried out by commercial software tools based on 
numerical methods. 

During the generic trade-off study, early in the design process, that is, when the concept 
of the spacecraft is still open, it is very convenient to assess the mean temperature of the 
satellite and its variation by means of analytical tools. This method provides valuable 
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information about the feasibility of the spacecraft from the thermal point of view as well 
as the possibility of easily carrying out parametric analyses to study the influence of the 
different variables involved in the problem. 

The purpose of the present study is to analyse the thermal behaviour of a small compact 
satellite spinning in a low orbit round a solar system planet (the Earth, say) as a function 
of time, using analytical or semi-analytical tools. In orbit, the spacecraft is subjected to 
periodic heat loads and it is expected to reach a thermally periodic state. Apart from 
obtaining the temperature oscillations once reached this periodic state, this analytical 
approach allows us to study the convergence of the satellite temperature to the periodic 
behaviour starting from different initial conditions, as it can occur from launch conditions 
or orbital manoeuvres. 

The results are obtained by applying the energy balance equation to the satellite, taking 
into account both the external heat loads, i.e. the solar irradiation, the planet albedo (solar 
irradiation reflected on the planet) and the infrared (IR) radiation from the planet, as well 
as the internal dissipation. 

We begin with the formulation of the energy balance equation as a non-dimensional 
ordinary differential equation (ODE), which is nonlinear and non- autonomous. Then, we 
consider an autonomous ODE that plays the role of an average of the actual equation. 
That autonomous equation is exactly soluble and guides our intuition of the behaviour of 
the actual equation and helps in its analysis. This analysis is carried out with qualitative, 
perturbation and numerical methods. Finally, we present our conclusions, regarding the 
design of satellites. 

2 Nonlinear ODE for the satellite temperature 

Let us assume that the satellite is approximately isothermal, so we characterize it by 
its temperature T. The heat input consists of an external flow coming from the space 
environment, in particular, the solar irradiation, the planetary albedo and the planetary 
IR radiation, and the internal heat due to the equipment dissipation. The solar irradiation 
and the planetary albedo have periodic variation, according to the position of the satellite 
on its orbit. On the contrary, the planetary IR radiation and the internal heat operate at 
constant rate. 

Let us assume that the solar irradiation changes from constant to vanishing (when 
the satellite is in the planet's shade), whereas the albedo depends on the angle between 
the sunlight and the vertical to the satellite (and is much smaller in absolute value). 
Furthermore, we assume that the albedo vanishes for one fraction of a period whereas the 
satellite's night (when it is in the planet's shade) lasts for another (smaller) fraction. We 
denote the heat rate due to solar irradiation by Qs, and the maximum albedo heat rate by 
Qa (at the satellite's noon). In addition, let the constant rate source have power Qc (the 
sum of the equipment dissipation and planetary IR radiation). Then, the energy balance 
equation that yields the temperature T is 

C f{t) = Qs fsii^t) + Qa faii^t) +Qc- Aea T{tf- 
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fs{x) = 1, < X < Xi or 1 — xi < X < 1; /s(x) = 0, Xi < x < 1 — Xi; 

fa{x) = COs(27rx), < X < X2 or 1 — X2 < X < 1; /a(x) =0, X2 < X < 1 — X2; 

fsA^) = fsA^ - 1), a; > 1- 

Here C is the satelhte's thermal capacity, u the orbital frequency, A is the satellite's surface 
area, e its emissivity, and cr the Stefan-Boltzmann constant. The values of xi and X2 are 
smaller than one half and they determine the fractions of the period with vanishing albedo 
or sunshine, respectively. In principle, X2 = 1/4 and xi > X2 (we assume that the albedo 
vanishes for half a period whereas the satellite's night is shorter: we take it as one fifth of 
the period, that is, Xi = 2/5, in an example below). 

We can write this equation in non-dimensional form by defining a = Ata / {Cv), kg = 
a}f^ Qs/ {Cv)^ ka = a^^^ Qa/iCu), kc = a^^^ Qc/ (Cu), and non-dimensional temperature 
variable 9 = a^/^ T and time variable ut (which we still denote t for notational simplicity). 
Then, 

m = kc + ksfsit) + kafait)-eit)\ (1) 

Unfortunately, this nonlinear ODE cannot be reduced to a quadrature; but we can deduce 
its relevant properties, nevertheless. 

We note that, if we remove the oscillating terms, the resulting equation 

e{t) = k- e{tf (2) 

(where k = kc) is immediately reduced to a quadrature, which can be integrated analyt- 
ically. Moreover, we can deduce the qualitative behaviour of the solutions in a straight- 
forward way: there is one fixed point, ^eq = k^^^, and it is stable. It is the equilibrium 
temperature, at which the heat input and the radiation output balance one another. When 
the temperature is close to it, 9 = k^^'^ + A^^, we obtain the linear ODE 

A9{t) = -Ak^/^A9{t), 

with solution 

A^(t) = A^(0) expi-Ak^^H). (3) 

Therefore, A9 halves in a time At = (In 2/4) k~^^^ = 0.17/6'eq^; and the larger is 9cq, the 
shorter it takes to reach equilibrium. We have checked that the exact solution of Eq. ([2]) 
is well approximated by the linear equation solution 

We remark that Eq. ([2]) can be naturally connected with Eq. ([1]) if we assume that k 
is the average of kc + kg fs{t) + ka fait) over one period, namely, k = kc + 2xiks + ka/n 
(rather than k = kc). We explore the consequences of this connection in the next section. 

With k = kc + 2xiks + /Ca/^r, reverting to physical variables, we have 



T 

eq 



1/4 

Qc + 2Xi Qs + Qa/'K 

Aea 



(4) 



Note that it does not depend on the parameters C and associated to the time derivative 
of the temperature in the ODE in physical variables. 
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3 Solution of the nonlinear ODE for the temperature 
evolution 



In this section, we first present a qualitative analysis, which allows us to prove the existence 
of an attracting limit cycle. To obtain the transient behaviour and the properties of the 
limit cycle, we employ perturbation theory. Finally, we perform a numerical analysis of 
the ODE, to check the results of the preceding methods and to obtain concrete and more 
precise numerical results. 

3.1 Qualitative analysis 

Equation ([T]) is non-autonomous and therefore equivalent to an autonomous system of two 
ODE's, namely, the system formed by Eq. ([T]) and the trivial equation t = 1. The generic 
behaviour of an autonomous system of two first-order ODE's is to have "simple" attracting 
sets, namely, equilibrium points or limit cycles, because of the Poincare-Bendixson theorem 
[2[ [3], H]. Given the nature of our problem, we can also apply the theory of ODE's with 
periodic coefficients (see, e.g., Ref. [5]): we can reduce the dynamics in the (t, 6')-plane 
to the cylinder [0,1) x (0, oo). Furthermore, we can consider 6 as the radial coordinate 
and 27rt as the angular coordinate of a plane with the origin excluded (equivalent to the 
mentioned cylinder). Given that there cannot be fixed points (since the equation t = 1 
forbids it), the Poincare-Bendixson theorem implies that the only possible attractors are 
limit cycles, with the period imposed by the heat input. 

The Poincare-Bendixson theorem states that a curve solution of a two-dimensional 
autonomous ODE that has no singularities (e.g., fixed points) and is contained in a compact 
domain for alH > approaches a limit cycle. We use a consequence of this theorem: given 
two concentric closed curves limiting an annular region that is free of singularities and such 
that the vector of derivatives on the two curves points towards the inside of the region, 
there exists a limit cycle in this annulus. 

To prove the uniqueness of the limit cycle, we appeal to Dulac's criterion for an annular 
region: if the autonomous system 

X = P{x,y), y = Q{x,y) 

is such that the divergence dxP + dyQ has constant sign in an annular region, then this 
region contains at most one limit cycle. The proof of this criterion follows from Green's 
theorem [3]. 

It is easy to check that Eq. ([1]) is such that 6 < 6^[^ = k^^^ =^ ^ > and 6 > 6^s.x = 
{ks + ka + kc)^^^ =^ 6 < 0, for all t. Note that 6*111111 and 6'max correspond to the equilibrium 
temperatures with constant minimum or maximum heat input, respectively. Therefore, the 
trajectories in the plane with 6 as the radial coordinate and 27rt as the angular coordinate 
which begin inside the annulus defined by those two temperatures are confined in it, and 
there is (at least) one limit cycle (with an oscillation in 6 confined to take place within those 
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values). Furthermore, this hmit cycle is unique, because div{kc + ks fs(t)+ka fa(t)—9^, 1) = 

The trajectories inside the annulus are alternately increasing and decreasing, as the 
right-hand side of Eq. ([T]) changes sign. To be precise, if we define the periodic function 

with minimum and maximum values O^i^ and 6^s.x, respectively, the sign of 6{t) changes 
when the trajectories cross it. This is represented in Fig. [H 

Using 6 as the radial coordinate and 2Trt as the angular coordinate. Fig. [1] becomes the 
polar plot in Fig. O This polar plot is more useful to represent convergence to the limit 
cycle in a standard way. 

We are interested in two questions: (i) finding the features of the limit cycle and the 
rate of convergence to it; (ii) analysing how this limit cycle and the rate of convergence 
depend on the constants kc,ks,ka. These questions cannot be answered by a qualitative 
analysis, so we turn to other methods. 

However, let us note that some more qualitative information can be obtained from a 
comparison with the non-oscillating Eq. ([2]): if O^q = k^^'^ = {kc + ^k^/^ + ka/i^Y^^ is 
sensibly smaller than one, the convergence time At = ^eq~^ is large, so that it is consistent 
to consider Eq. ([2]) as an averaged equation. To be precise, Eq. ([2]) can be derived as an 
equation for the mean temperature in Eq. ([T]) by averaging it over one period. Therefore, 
the evolution is given by a long-time decay to ^cq and a short-time oscillation about that 
long-time behaviour (as in Fig. [1]). 

Thus, it seems convenient to consider Eq. ([1]) as a perturbation of Eq. ([2]), namely, to 
consider the time-dependent (oscillating) functions /a s as a perturbation. Then, we can 
employ standard perturbation methods [Hj. We do so in the next section. 

3.2 Perturbation theory 

As we have explained above, our intention is to take advantage of the simple solution of 
Eq. ([2]) by taking it as the basis of a perturbation scheme. Therefore, we write Eq. ([T]) as 

e{t) = k + ef{t)~9{t)\ (5) 

where k = kc+2xi kg + ka/n, f{t) = kc — k+kg fs{t) + ka fa{t), and e is a formal perturbation 
parameter to be set to one at the end. Note that f{t) is defined such that it has vanishing 
average over a period and, therefore, it represents the deviations about the mean value 
k. Furthermore, the solution of Eq. ([5]) must fulfill the initial condition 9{0) = 9in- We 
assume an expansion of 6{t) of the form 

oo 

e(t) = J]e"e„(t), 

n=0 
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temperature evolution 




Figure 1: Results of the numerical analysis of the ODE with = 0.13, ka = 0.007 and 
kc = 0.016, explained in the text: (left) convergence to the attractor; (right) detail of the 
first period. Note how the trajectories change from increasing to decreasing, and viceversa, 
as they cross the dashed line 6\i^{t) (the bottom of this line, at kj'^'^ = 0.36, is below the 
plot regions). 




Figure 2: Polar plot of the trajectories in Fig. [T] (right): the radius is 6* — 0.5 (subtracted to 
display well-separated lines) and the angle is 27rt. The limit cycle is the bold dashed line. 
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and substitute it into Eq. ([5]). Equating to zero the successive powers of e, we obtain the 
basis equation and an infinite set of hnear ODE's, namely, 

eoit) = k-eoit)\ (6) 

9,it) = fit)-i9o{t)%{t), (7) 
9,{t) = -69o{tfe,{tf - 49o{t)'92{t), (8) 

The initial condition for 9o(t) is 6*0(0) = 6'in, while the remaining equations fulfill 6'„(0) = 0. 

These ODE's can be solved subsequently, namely, we can solve first the equation for 
^o(^) and, then, we can solve the equation for 9i(t) and so onwards. Of course, the equation 
for ^o(^) is Eq. ([2]), which we have already solved. The following equations are just first- 
order linear inhomogeneous equations, which are soluble by quadratures [6J. However, their 
solution involves complicated integrals that cannot be made analytically. In particular, 
Eq. (iTl) has the solution 

ei{t) = J^l fir)I{T)dT, /(t) =exp 4^ 9oirfdr ; 

but I{t) has no analytic expression. 

To overcome this difficulty and obtain analytic expressions, we can use instead of the 
exact function 6^0 (t) the approximated function given by Eq. ([3]), namely, 9o(t) = 9cq + 
{9in — 9eq) exp(— 4/^^/"^^) . Thus, the integral in I{t) becomes a sum of integrals of exponential 
functions, which yield exponential functions again. In the large-t limit (when t ^ 9^^^), 
the exponentials decay and the value of I{t) further simplifies to I{t) = exp (4/c'^/*^t) . Then, 

9^{t) = [ f{T) exp [-4k^/\t -T)]dT= [ fit- t) exp [-Ak'^/^r] dr. (9) 
Jo Jo 

Furthermore, in the large-t limit and given that f{t) is periodic, we can extend the upper 
integration limit in the latter integral from t to oo. Thus, we have a periodic function, 
namely, an approximation to the limit cycle (past the transient regime). We can obtain a 
general expression of this periodic function through Fourier analysis. 

3.2.1 Fourier analysis 

Let us expand /(t) in Fourier modes as 



fit) = E 



m 

m=—oo 



Since f{t) is real and we assume that it is symmetric with respect to t = 0, the coefficients 
frri are real and f-m = fm] in addition, /o = fit) dt = 0. The other Fourier coefficients 
are given by 

fm= I /(t)e-2-™*rft. (10) 
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It is straightforward to solve for 6i{t): 



^ ^ fm e^™* _ 9 f 4^^^^ cos(27rmt) + 27rm sin(27rmt) 

m,=— oo m=l 

If we substitute the form of f (t) = kc — k + kg fs(t) + ka fait) into Eq. ffTOj) . we obtain 



_ sin(27rmxi) kg 

Jm Kg "To 

vrm 2 
when m 7^ 0. With X2 = 1/4, 



sin(27r(m + l)x2) ^ sin(27r(m — l)x2) 



7r(m + 1) 7r(m — 1) 



sin(27rmxi) /c^ / tt 

Jm = kg 7 — ^ r COS m — 

nm nym^ — 1) V 2 

On the other hand, independently of the form of fm, if the convergence time At = ^eq^^ 
is large, we have that Ak^/^ <C 27r, so we can neglect in Eq. ( fTTl) the terms with k and write 
the crude approximation 

^ / X „ sin(27rmt) 

m=l 

This approximation is equivalent to neglecting in the right-hand side of Eq. ([7j) the second 
term with respect to the first one, that is to say, it implies that 9i{t) follows f{t) with no 
delay. We can use it to provide an estimation of the amplitude of the oscillations, with 
the following simple method. This estimation can be useful in the conceptual design of the 
satellite. 

The largest variations of the slope of 9{t) and of its approximation 9i(t) take place at 
the discontinuities of /(t), namely, at t = xi, 1 — xi (in the first period). f{t) changes 
sign at those times and the periodic function 9i{t) has its maximum and minimum there. 
In between, we can take f{t) constant and, therefore, 9i{t) linear. Hence, we compute the 
slopes at two convenient points, namely, t = 0, 1/2, where 9i{t) = 0, according to Eq. f[T^ . 
We use this information to calculate the maximum and minimum of 9i{t). The slope at 
t = is /(O) = (1 — 2xi)ks + (1 — l/7T)ka, so the maximum is 

^Imax = /(O) = Xi [(1 - 2xi)A;, + (1 - 1/Tr)ka] . (13) 

Then, with kg = 0.13, ka = 0.007 and kc = 0.016 (values to be justified in Sect. 13.31) . 
^imax(a;i = 0.4) = 0.4 X 0.031 = 0.012. The slope at t = 1/2 is /(1/2) = -2xiA;, - ka/rr, 
so the maximum is 

^imax = -(1/2 - Xi) /(1/2) = (1/2 - Xi) [2x,kg + ka/n] . (14) 

In particular, 6'ijnax(3;i = 0.4) = (0.5 — 0.4) x 0.106 = 0.0106. Both values agree sufficiently. 
Since 9i(t) in Eq. (fT2l) is an odd function, the minimum is at t = 1 — Xi = 0.6 and its value 
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is the negative of the maximum. A graphical comparison of the present approximation 
with the actual limit cycle is displayed in Fig. [31 

Note that Eqs. ( fT3ll and (fT4l) coincide if fc^ = 0. In fact, the preceding method relies on 
the particular form of f{t), namely, it is close to a step function. This is the reason why 
the method yields a good result. Of course, while still using the approximation leading to 
Eq. f[T2|) . the preceding method is improved by the formula 



which demands a little more work to yield the maximum value of 6i. Better approximations 
to 6*1 are possible by considering k in Eq. (fTTl) . that is to say, by considering the displacement 
of Oi{t) [Oi{0) 7^ 0]. In a different sense, the approximation is also improved by carrying 
on the perturbation scheme to the second order. 

3.2.2 Second order approximation 

Eq. ([8]) is similar to Eq. ([7]), if we consider Oq = 9^^ and (the square of) the first order 
solution as the forcing term. Therefore, a similar reasoning leads us to the existence of a 
periodic solution for 62{t) in the large-t limit, which can be obtained by Fourier analysis. 
The result is 



We can see that the perturbation series is a power series in the Fourier coefficients of 
the forcing term and its convergence depends on the magnitude of this term. Given that 
this term is proportional to kg and ka, the convergence of the series is improved when these 
constants are small. This condition is related to the one that we found for Eq. ([2]) to hold 
as an average equation, namely, that k = kc + 2x\ kg + ka/7r is small. Thus, it is consistent 
with Eq. as the basis of the perturbation scheme. 

3.3 Numerical solutions 

We use a standard numerical method to integrate Eq. ([1]), for given values of ks,ka, kc- Let 
us see what values are adequate for a small satellite that is orbiting round the Earth and 
which can be considered isothermal. 

We can take the satellite to be a cube of 0.5 m side and mass of 50 kg. Let us 
assume that it is mostly covered by solar cells. Then, suitable values for the absorptivity 
and emissivity are a = 0.8 and e = 0.7, respectively. We assume the satellite's average 
specific thermal capacity to be 0.2 W h kg~^ (corresponding to a composition of 
alluminium with some plastic). An adequate value for the orbital frequency is (1.5 h)~^. 



Therefore, Cu = 50kg 0.2 Whkg"^ K-i(1.5 h)-i = 6.67 W/K and a = Aea/{Cv) = 
1.5m2 0.7 5.710-*^Wm-2K-V(6.67W/K) = S.ggiO^^K-l 





p= 



(15) 



9 



linear approx. 




t 



Figure 3: Linear approximation to the first-order perturbative limit cycle 6'eq + Oi(t), such 
that it is linear by pieces (dashed line), compared to the real limit cycle (full line). 

The solar irradiation heat input is the product of the solar constant, the absorptivity 
and the projected area (a quarter of the real area); namely, Qs = 1370 W m~^ 0.8 [6(0.5 m)^/4] = 
411 W. Therefore, fc, = a^/^Qs/{Cu) = 0.128. The albedo is very variable, and to calculate 
Qa, it is necessary, in addition, to consider the visibility factor. We take an average value 
of Qa = 0.05 Qs so ka = 0.0065. 

The constant heat input is the sum of the equipment dissipation plus the planetary IR 
irradiation input. The equipment dissipation power is due to the transformation of the 
incoming irradiation power. If we assume that the solar cells cover 80% of the satellite's 
surface and their efficiency is about 10%, and we take into account that the solar irradiation 
only holds for a fraction 2xi = 4/5 of the period, we deduce that the power dissipated is 
about 30 W. The planetary IR irradiation input can be estimated by applying the energy 
balance to the planet (Earth). This balance yields an IR radiation flux of about 7% of 
the solar constant. We can take the absorptivity in the IR to be similar to the emissivity 
e = 0.7. Hence, the planetary IR irradiation is about 23 W and Q^ = 53 W. 

Summarizing, we can take = 0.13, ka = 0.0070 and k^ = 0.016. Furthermore, we 
consider seven values of the initial condition T(0) (t = is the noon) uniformly distributed 
in the interval 273 K < T < 300 K, corresponding to nondimensional 0.57 < 9 < 0.63. 
These values are included between ^min = kc^'^ = 0.36 and 6'max = (kg + ka + kc)^^^ = 
0.625. The results are plotted in Fig. [H for the time interval < t < 5, which is sufficient 
to show convergence to the limit cycle. 

Indeed, in our example, we are in a situation in which 6'eq = {k^ + 4/^^/5 + ka/i^Y^'^ = 
0.591 is sufficiently small to apply the averaged equation ([2]). The convergence time is 
At = 9eq~^ = 4.8. Fig. [1] clearly shows that the evolution is given by an oscillation with 
unit period superimposed to a slower and approximately exponential convergence in a time 
At. 

To measure precisely the attractor characteristics, we choose the initial condition 6{0) = 
6eq = 0.591 and we let the numerical integration up to t = 10. The resulting solution has 
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0(10) = 0.592, which is shghtly different from ^eq (although the difference is inappreciable 
with two-digit precision). Of course, the solution with 9{0) = 0.592 is the limit cycle. 
It has local minimum and maximum amplitudes 0.58 and 0.60 at t = 2/5 and t = 3/5, 
respectively. Those extreme values are quite close, considering the interval defined by 
6*111111 = 0.36 and 6'max = 0.625 (see Fig. [1]). The limit cycle is plotted in Fig. [3l 

4 Discussion 

Regarding the problem of heat transfer in a satellite, the main conclusions that we can draw 
from our analysis are the following. First of all, as regards the temperature stability, it 
holds in the sense of convergence to the limit cycle behaviour, guaranteed by the application 
of the Poincare-Bendixson theorem and Dulac's criterion for an annular region. Moreover, 
the convergence is exponential. 

The limit cycle behaviour has a time dependence related to the periodic heat input, 
namely, approximately similar to its integral, but somewhat displaced. To be precise, the 
limit cycle consists, in one period, of a temperature-growing phase when the solar heat input 
is on (as is the albedo), and a relaxation phase in which the temperature falls to adjust to 
the constant input Qc (although this phase ends long before the temperature approaches 
the corresponding equilibrium value). The total amplitude of the temperature oscillation 
is relatively small in the example that we have studied. This example, with realistic values 
of the non-dimensional constants kg, ka, k^ for a low-orbit small satellite (reasonably small 
values, in particular), shows that the attractor is approached exponentially but in a time 
reasonably larger than the period. This time is simply given in terms of kg, ka, kc by 
At = k~^^^, where k = kc + 2xiks + ka/Ti- Moreover, the mean temperature in the limit 
cycle is 9eq = k^^^. 

It is useful to compare the mean temperature ^eq with the oscillation about it. The 
former depends on fc, namely, the average of kc + kg f sit) + kafait), whereas the latter 
depends on f(t), which is the oscillation of ksfsit) + kafait) {kc is a constant). Thus, the 
temperature oscillation is independent of kc (the constant heat input). This constant is 
bound to be quite smaller than kg, for physical reasons, since the equipment dissipation 
power is due to the transformation of the incoming solar irradiation power, and the IR 
irradiation input is a fraction of the solar input. On the other hand, ka <C kg, and we 
assume that Xi > 1/4. In conclusion, the main contributions to k and, therefore, to ^eq 
are due to ks. If the temperature oscillation is estimated by Eqs. ( fT3|) or ( IT4l) . the main 
contribution to it also seems to be due to kg] except for the factor 1 — 2xi, which can be 
small. Let us first note that Eq. ([T^ yields a larger value than Eq. ([H]) if xi > l/(27r), 
which always holds. Focusing on Eq. (fT3|l . we note that the coefficient 1 — 2xi of ks is 
smaller than the coefficient of ka- When xi approaches 1/2, we only have the fca-term: then 
the satellite is always under the solar irradiation, which does not oscillate; in consequence, 
the temperature oscillation is only due to the albedo and, therefore, is depressed. 
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Let us express the temperature oscillation given by Eq. (fT3!) in physical variables: 



Contrary to the expression of Teq in Eq. (jl]), this expression depends on C and v but it does 
not depend on A, e and a. We expect both expressions to be very useful for conceptual 
thermal design. 

Of course, in the numerical example in Sect. 13.31 we have not considered every orbital 
circumstance and we have just intended to find a set of sensible values for the constants. 
More information on satellite design and the space thermal environment can be found in 
the literature [T]. It is worthwhile to discuss briefly here possible changes of the values 
of ks and kc- While the ratio k^/kg ~ 0.1 is adequate for a small satellite, the absolute 
values of ks-,kc can be amply changed; for example, by increasing the period (which we 
have taken as 1.5 h in our case). It is easy to see that both k^ and k^. are proportional to 
the 4/3rd-power of the period. Thus, if we increase the period by a factor of ten, say, the 
constants increase by a factor of 10^/^, so that At = (fee + 2xiks + ka/T^)~^^^ decreases by a 
factor of ten, becoming smaller than one (the period). In this situation, the convergence to 
the hmit cycle is very fast. In fact, it is appropriate to consider the opposite approximation 
to the temperature evolution, namely, during the decay to the limit cycle, the temperature 
only changes slightly along this limit cycle. 

Further analysis could be based on a more refined model than the isothermal model that 
we have used. In particular, a two-node model has been studied by Guerra, Perez-Grande 
and Sanz [7|. The two ODE's for the two-node model are equivalent to an autonomous 
system of three ODE's. As is well known, such system can have chaotic behaviour [4]. 
However, we expect that the situations that are well approximated by the one-node model 
that we have analysed do not exhibit chaotic features and, therefore, can be conveniently 
studied, and provide, in addition, a basis for this very simple model. 
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